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Abstract 

The problem of gravity propagation has been subject of discussion for 
quite a long time: Newton, Laplace and, in relatively more modern times, 
l Eddington pointed out that, if gravity propagated with finite velocity, 

O planets motion around the sun would become unstable due to a torque 

O"* originating from time lag of the gravitational interactions, 

j.^ Such an odd behavior can be found also in electromagnetism, when 

W) one computes the propagation of the electric fields generated by a set 

of uniformly moving charges. As a matter of fact the Lienard-Weichert 
I retarded potential leads to a formula indistinguishable from the one ob- 

tained assuming that the electric field propagates with infinite velocity. 
Feyman explanation for this apparent paradox was based on the fact that 
| uniform motions last indefinitely. 

^\ To verify such an explanation, we performed an experiment to measure 

04 the time/space evolution of the electric field generated by an uniformely 

I moving electron beam. The results we obtain on such a finite lifetime 

kinematical state seem compatible with an electric field rigidly carried by 
£SJ the beam itself. 

1 Introduction 



X 



In Space, Time and Gravitation Eddington discusses [I] the problem of 
gravity propagation. He remarks that if gravity propagated with finite 
velocity the motion of the planets around the Sun would become unstable, 
due to a torque acting on the planets. The problem was already known 
to Newton and was examined by Laplace [2] who calculated a lower limit 
for the gravity propagation velocity finding a value much larger than the 
speed of light. 

However, at the time of Eddington's writing, General Relativity had 
been just formulated, with gravitational waves travelling with the speed 
of light as a possible solution. Eddington noted that a similar problem 
existed in electromagnetism, and since electromagnetic waves in vacuum 
do travel with the speed of light, he concluded that in General Relativity 
gravity also propagates with the speed of light. 
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We remark that an intriguing behavior of electromagnetism occurs 
computing the field of an electric charge moving with constant velocity. 
One finds that in such a case the electric field at a given point P(x, y, z, t) 
evaluated with the Lienard-Wiechert potentials is identical to that cal- 
culated by assuming that the Coulomb field travels with infinite velocity. 
This is explained by Feynman .3; assuming that the uniform motion occurs 
indefinitely. 

However, it is important to consider that the Lienard-Wiechert (LW.) 
retarded potentials represents just one particular solution of the Maxwell 
equations. The general solution is represented by the sum of the retarded 
as well the anticipated potentials, the last one usually dismissed because in 
conflict with the causality principle. If we consider the general solution, we 
find that the electric field may, in addition to the well-known propagation 
with the speed of light, also propagate instantaneously. The only way to 
shed some light on this problem is by means of an experiment. 

To verify if the Feynman interpretation of the L. W. potentials holds in 
case of a charge moving with constant velocity for a finite time, we have 
performed an experiment to measure the time evolution of the electric 
field produced by an electron beam in our laboratory; such kinematic 
state has obviously a finite lifetime. 

It is well known that a sizeable number of instrumentation devices 
(e.g. beam position monitors) are based on effects produced by electric 
fields carried by particle beams. The effects and the propagation of such 
fields, however, have never been studied in details: the main point ex- 
ploited by these devices is that the field effects are contemporary to the 
particles passage and that the signal size obtained, for instance, on a pair 
of striplines inside a vacuum pipe yield a (transverse) position measure- 
ment of the beam itself. The experimental situation for those devices is 
quite complicate to understand, as all the fields are inside a conductor and 
the transverse distanced exploited are always small. We, on the contrary, 
tried to carry out our experiment in a clean environment: the electron 
beams used were propagating in a vacuum like environment. We covered 
a wide range of transverse distances w.r.t. the beam line (up to 55 cm). 
Such range leads to explore time and space domains for the emission of 
the detected field far outside the physical region. 

We have found that, in this case, the measurements are compatible 
with an instantaneous propagation of the field. We believe that this in- 
triguing result needs a theoretical explanation in addition to the Feynman 
conjecture or the naive hypothesis of instataneous propagation. 

2 Theoretical Considerations 

The electric field at r(x, y, z) from a charge e travelling with constant 
velocity v, at a time t can be written, using the Lienard-Wiechert retarded 
potentials as [U |B] : 
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where 

R{t') =r-vt' (2) 

is the distance between the moving charge and the space point where one 
measures the field at time t, and 

f = t-W (3) 
c 

The field from a steadily moving charge can also be written (as easily 
deducible from eqn. [I] in case of constant velocity) [3J [4j [5j [6] as 

_e_ *W 1-^ 

U 47r eo ii(*)3 (1 _£ sin 2 ( (t)) ! lj 

where J2(t) is the vector joining the charge position and the point at which 
we evaluate the e.m. field at time t (Eqs. 38.8 and 38.9 of[D)[]]and 9 if) 
is the angle between v and R(t) 

A pictorial view of the above mentioned quantities can be seen in Figjl] 




Figure 1: A pictorial view of various quantities mentioned in eqns [T] and [4] 

If we indicate with y the generic transverse coordinate, using eqn[l]we 
can compute the maximum transverse electric field w.r.t. the direction of 
motion, given by (7 = 1/ y/1 — v 2 /c 2 ) : 

47re y J 



1 In Landau's words: the distance R(t) at precisely the moment of observation (see 

pag. 162 in g]? ). 
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a value obtained when the charge is at a distance ■yy at a time 

t' = t-™ ( 6 ) 
c 




Figure 2: The electric field from eqn.[T] normalized to its maximum value, 
E y (R)/E max , generated by 500MeV electrons as a function of z' (or t', lower 
abscissa scales), expected at (z = Ocm, y — 30 cm). z' and t' are defined in 
ecms.[2] and [3j The horizontal scale of the upper graph (a) is such to include 
the point where E y (R) = E max ; the lower graph (b) is a close-up of the region 
z G [—10, 0] m typical of our experiment (note the different vertical scales). 



Fig. [2] shows the field, normalized to E max , generated by relativistic 
electrons (E — 500 MeV) moving along the z axis, at a transverse distance 
y — 30 cm. We observe that the maximum value of the field appears 
to be generated when the charges are in an unphysical region, namely 
z = —300 m. Conversely, the electric field in the region of our experiment 
(|z| < 10 m) is many orders of magnitude smaller. 
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Figure 3: A sketch of our experimental apparatus (side view). Sensors A2, A4, 
ortogonal to the figure plane are not shown. 



In our experiment we measure the electric field generated by the electron 
beam produced at the DA5>NE Beam Test Facility (BTF) [7J, a beam 
line built and operated at the Frascati National Laboratory to produce 
a well-defined number of electrons (or positrons) with energies between 
50 and 800 MeV. At maximum intensity the facility yields, at a 50 Hz 
repetition rate, 10 nsec long beams with a total charge up to several 
hundreds pCoulomb. The electron beam is delivered to the 7 m long ex- 
perimental hall in a beam pipe of about 10 cm diameter, closed by a 40 /im 
Kapton window. Test were carried out shielding the exit window with a 
thin copper layer, but we did not observe any change in the experimental 
situation. At the end of the hall a lead beam dump absorbs the beam 
particles. In our measurements we used 500 MeV beams of 0.5 ■£■ 5.0 x 10 s 
electrons/pulse. (7 ~ 10 3 ). 

A schematic view of the experimental setup is shown in Fig. [3] At 
the beam pipe exit flange, electrons go through a fast toroidal trans- 
former measuring total charge and providing redundancy on our LINAC- 
RF based trigger. 

To measure the electric field we used as sensors 14.5 cm long, 0.5 cm 
diameter Copper round bars, connected to a FLASH-ADC by means of a 
fast, terminated coax cable. 

To record the sensors waweforms we used a Switched Capacitor Array 
(SCA) circuit (CAEN mod V1472) able to sample the input signal at 5 
GSamples/sec. In addition to the sensors output, the SCA stored also the 
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LINAC-RF trigger and the toroid pulse. 

The Coulomb field acts on our sensor quasi-free electrons, generating 
a current. An example of the recorded signals is shown in Fig|4] The 
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Figure 4: Typical A5 (left) and A6 (right) sensor responses. Lower plots show 
in detail the granularity of our time measurements. (200 psec./bin) 

pulse shape, unlike the current intensity, depends on the inductance, ca- 
pacitance and resistance (L, C, R) of the detectors. 

The sensor response V(t) for a step excitation Vb can be written as: 

V(t) = V -e-£ t sin(ujt). (7) 

The voltage difference between the bar ends for the maximum value of 
the Coulomb field, obtained suitably modifying eqn. [5] for a finite longi- 



tudinal extent of the charge distribution, (in our case, the electron beam 
is « 3 m. long) is: 

r ,t A (y+ 14.5 cm\ 

Vmoa = ??" Zn , (8) 

27re D \ V J 

where A is the charge per unit lenght of the incoming beam and r\ is the 
sensor calibration constant. In the electric field calculations, the image 
charges appearing on the flange as the beam exits the pipe have also been 
included. However, as their effect decreases rapidly with the distance from 
the flange, it is completely negligible in our experiment (distance^ lm). 
The sensor calibration has been carried out using a known field generated 
by a parallel plate capacitor. We find experimentally r\ — 7.5 x 10~ 2 ± 3%, 
however due to various systematic effects we believe our calibration to be 
good to « 20%, in absolute terms. 

Assuming that the L.W. formula (eqn [8| holds (which should apply 
only if the uniform charge motion would last indefinitely and the charges 
generating the field would not be shielded by conductors) we expect, in 
our typical beam operating conditions, pulse height of the order of 10 mV 
out of our sensors. In the more realistic hypothesis that the L.W. formula 
should be corrected to take into account the beam pipe shielding and the 
finite lifetime for the charges uniform motion, as it is in our experiment, 
the expected amplitude, cfr. fig|2j would be of the order of few nanoVolt 
and hence unmeasurable. 

We used six sensors: four of them, A1,A2,A3 and A4 in the following, 
located at a (longitudinal) distance of 92 cm from the beam exit flange 
(cfr. fig [5]), in a cross configuration, each at a transverse distance of 5 cm 
from the beam line. The main purpose of these four sensors is to provide 
reference for the other two detectors A5 and A6 located through out the 
measurements at various longitudinal and transverse coordinates along 
the beam trajectory. 



4 Measurements and data base 

Electron beams were delivered by BTF operators at a rate of few Hertz; 
data were collected in different runs, identified by given longitudinal and 
transverse position of the movable detectors. (A5 and A6) 

We collected a total of eighteen runs, spanning six transverse positions 
and three longitudinal positions of A5 and A6 for a total of w 15,000 trig- 
gers. Through out the data taking, the references sensors (A1,A2,A3,A4) 
were left at the same location in order to extract a timing and amplitude 
reference. As mentioned before, we collected data with the movable sensor 
at 172 cm, 329.5 cm and 552.5 cm longitudinal distance from the beam 
exit flange. For each of the longitudinal positions we collected data on six 
transverse positions: 3, 5, 10, 20, 40 and 55 cm from the nominal beam 
line. For each run, A5 and A6 were positioned symmetrically with respect 
to the nominal beam line; spatial precision in the sensor positioning was 
« few mm in the longitudinal coordinate m 1 mm. in the transverse. 

We define: 

Sn = Vma * X (9) 

1 * elec, 
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where V ma x is the peak signal recorded by the SCA and N e i ec . is the total 
number of electrons in the beam, as measured by the fast toroid. The 
10 s factor in eqn[9] takes into account the typical beam charge. As an 
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Figure 6: Normalized amplitude Sn (see text for details) for sensor Al in two 
different runs: the above mentioned sensor was in the same operating conditions 
through out the different runs, taken minutes apart. 
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example, the two plots in Fig. [6]show the Sn values for the reference sensor 
Al, for two different runs. One issue common to all our measurements 
stands out clearly: in the same experimental conditions (sensor position, 
trigger timing, cable lengths, DAQ settings) the two distributions are 
different, We attribute this difference to less than perfect reliability in 
the beam delivery conditions (launch angles, total beam length, charge 
distribution in the beam pulse length, stray magnetic fields, etc.), over 
which we had little control. 

Since our four reference sensors must provide normalization for the 
measurements taken by A5 and A6, our analysis proceeds as follows. We 
obtain, on a run-by-run basis, the average amplitude of the four reference 
sensors, either evaluating the medians, or fitting the distributions with 
gaussian functions and taking the mean values. 

Next we make the assumption that in any given series of runs under 
study, different only for the position of the movable sensors A5 and A6, 
each reference sensor, that is never moved, must always yield the same 
amplitude. Since Fig. [6] shows that this is not the case, we need to allow 
for some (uncontrolled) effect due to variation of beam parameters. We 
do so by enlarging the errors on the reference amplitudes, originating from 
the gaussian fits, by a rescaling factor. This factor is chosen by requesting 
that the reduced \ 2 of the series be consistent with the hypothesis that 
all amplitude measurements for any given sensor in the series have one 
common value. 

Once obtained the error rescaling factor, for a given sensor and run 
series, we proceed to analyze the movable detectors by enlarging by the 
same factor their own uncertainties. The amount of rescaling needed is of 
the order of 10; the overall relative error on the run-averaged pulse height 
is typically w 10%. 

4.1 Amplitude as a function of transverse dis- 
tance 

As mentioned in the previous section, at each longitudinal position we 
collected data at six different transverse positions for A5 and A6. The 
requirements placed on data were: a lower cut on the beam charge (N e > 
0.5 x 10 s ) and upper cut on the baseline noise on the six detectors (noise 
< 0.5 mV , where the typical r.m.s. noise was 0.15 mV). We select in this 
way roughly 70% of the recorded triggers. 

Fig. [7] [5] and [9] show normalized amplitudes Sn versus transverse dis- 
tance obtained for the three different longitudinal positions. The displayed 
results are completely consistent with eqnjHJ which gives the voltage value 
expected in case of a charge indefinitely moving with constant speed. 



The results shown in fig |7|8| and [9| were obtained without any normal- 
ization between measurements and LW. theory. We stress again that the 
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Figure 7: Normalized amplitudes Sn for sensors Al(upper) and A3(lower) 
for za5,A6 =172.0 cm. The continuous line at 8.98±0.54 mV (upper) and 
10.25±0.59 mV (lower) indicate the weighted average of our measurements. 
The six measurements plotted refers to the transverse positions of sensor A5, 
A6. No dependence is expected (see text), as sensor Al, A3 were kept in the 
same operating conditions and at the same locations. The dashed line indicates 
the nominal normalized value V^ax — 9-78 mV of eqnj8]for y=5 cm. 
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Figure 8: The same plot as in figj7] at za5,ag= 329.5 cm. The continuous lines 
at 8.80 ±0.63 mV (upper) and 10.47±0.86 (lower) indicate the weighted average 
of our measurements. 
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Figure 9: The same plot as in figj7]at ZA5,;16 — 552.5 cm. The continuous line at 
9.67 ±0.54 mV (upper) and 10.75±0.74 (lower) indicate the weighted average 
of our measurements. 
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amplitude we measure is many orders of magnitude higher than the one 
would expect from the unshielded beam charge. Were we sensitive only to 
fields generated by the electron beam once they exited the beam pipe, our 
pulse height would have been, as mentioned in the previous paragraph, in 
the few nanoVolt range and then undetectable. 

In figg. |10| 1 1 1 1 and [12] we show the amplitude ratios between sensors 
A5 and Al (A6 and A3) as a funtion of transverse distance from the beam 
line. Also in this case, data are completely consistent with the logarithmic 
behavior of eqn [8] 




10 20 30 40 50 60 



cm 

Figure 10: Upper graph: the points show the ratio ^/" ax |^ at z j 4 5 j 4 6 =172.0 cm 
versus the transverse distance. Lower graph: y""" [^3] . The continuous lines 
represent eqn. [8] for the depicted ratios. The two reduced % 2 are respctively 
1.82 and 1.06. No fit has been performed on the data: the reduced x 2 has been 
evaluated from eqn [8] and the experimental data. 

4.2 Timing measurements 

Our 200 psec/chn SCA provides timing for detector outputs, so that it is 
possible to detect both longitudinal and transvese position-time correla- 
tions. As a reminder we stress again that, in the hypothesis of stationary 
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constant speed motion, no time difference is expected as a funtion of trans- 
verse distance, while different longitudinal positions should exhibit delays 
consistent with particles travelling at 7 w 1000. 

Also in this case, we will have to rescale the errors yielded by standard 
procedures extracting central values from quasi gaussian distributions; in 
this case we impose that, by symmetry, the time difference between A5 
and A6 be independent from the transverse distance between detector 
and nominal beam line. We then duplicate the procedure described in 
the previous section requiring that time difference between A5 and A6 for 
each run all come from a common value. 

The upper graph of fig [13] shows the time difference relative to 172 cm 
longitudinal distance data. The amount of rescaling, in this case, is about 
a factor of 10 and the overall resolution on time difference measurements 
is of the order of 50 psec. 

The data show no time dependence of the sensor signal on transverse 
distance: the reduced x 2 for the hypotesis of a constant delay as function 
of y is always below 2 at each longitudinal positions. Furthermore, would 
one add a linear term depending on tranverse distance for the sensor 
time delay, the inverse velocity obtained would have a value smaller than 
3 x 10" 9 ^ at 95% confidence level. 

We summarize the time distance correlations in the Table [l] where the 
data obtained at the three different longitudinal positions are shown. 



Table 1: Timing measurements. The expected differences are calculated for 500 
MeV electrons. 



longitudinal distances 
between two sensors [cm] 


expected experimental experimental 
[ns] a5 [ns] a6 [ns] 


(552.5-329.5) 223.0±1.5 
(552.5-172.0) 380.5±1.5 
(329.5-172.0) 157.5 ±15. 


7.43 ±0.05 7.28 ± 0.02 7.52 ± 0.04 
12.68 ± 0.05 12.62 ±0.04 12.84 ±0.05 
5.19 ±0.05 5.21 ± 0.03 5. 17 ±0.04 



5 Discussion 

The highly relativistic beam we have used does not allow to distinguish 
on a timing basis between the effects of the electron beam and any elec- 
tromagnetic disturbance generated e.g. at the beam pipe exit windows. 
In principle such a discrimination might be possible using low energy par- 
ticles, with /3 <C 1. We note, however: 

1. Any e.m. radiation could only originate from a point inside the 
beam pipe, as there is no material in the beam line and no relevant 
electric/magnetic field is present in the experimental hall after the 
beam exit flange. This implies that, moving away from the exit 
window, e.m. disturbances will eventually cause smaller and smaller 
responses from the detectors, which is not shown by figg[l0] [TTJand 

mi 
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Figure 13: Top graph: time difference between A5 and A6 versus the transverse 
distance. Errors have been rescaled according to the procedure described in the 
text. Middle and bottom plots: time difference between the movable sensors 
A5 and A6 (respectively) and one of the fixed sensors, Al. ZA5.A6 = 172 cm. 
The line at 1.549±0.036 nsec. (middle) and 1.500±0.036 nsec (lower) indicate the 
weighted average of our measurements, once the reduced % 2 is rescaled according 
to the procedere described in the text (par. 3.2). 
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2. It is possible to control the horizontal beam launch angle in the ex- 
perimental hall varying the current setting of the last dipole magnet. 
One can then predict the amplitude ratio of two detectors located 
right and left of the beam line. Special runs were taken to this 
purpose and the results are completely consistent with the expected 
variation of the horizontal beam displacement. 

3. At the moment it is not possible to exclude that low frequency 
Cherenkov and/or transition radiation effects [S] might contribute 
to the detected signal. We plan, in the near future, to carry out the 
experiment with an evacuated non conductive beam pipe installed at 
the end of the beam delivery flange and moving the kapton window 
downstream with respect to our sensors, in order to reduce the two 
above mentioned effects pratically to zero. 

With reference to table [T] we notice that the longitudinal time differ- 
ences are completely consistent with the hypotesis of a beam travelling 
along the z axis with a Lorentz factor 7 w 1000. 

Such an occurence agrees with the Lienard-Weichert model. Retarded 
potentials, however, predict that most of the virtual photons [§] responsible 
for the field detected at coordinates z and y be emitted several hundred 
meters before the sensor positions and at different times according to the 
detectors transverse distances. Conversely, assuming that such virtual 
photons are emitted in a physically meaningful region (between the beam 
exit window and our detectors), the amplitude response of the sensors 
should be several order of magnitude smaller than what is being measured 
(cfr. Figj5|. Our result, obtained with a well definite set of bounday 
conditions (longitudinal and transverse distance between beam line and 
sensors, details of the beam delivery to the experimental hall etc.) matches 
precisely (within the experimental uncertaities) the expected value of the 
maximum field calculated according to L.W. theory, that is also the value 
calculated with eqn|4] when the beam is at the minimum distance from 
the sensor. 

We again point out that the consistency of our measurement and eqn 
[8] has been obtained without any kind of normalization. 

6 Conclusions 

Assuming that the electric field of the electron beams we used would act 
on our sensor only after the beam itself has exited the beam pipe, the L.W. 
model would predict sensors responses orders of magnitudes smaller than 
what we measure. The Feynman interpretation of the Lienard-Weichert 
formula for uniformly moving charges does not show consistency with 
our experimental data. Even if the steady state charge motion in our 
experiment lasted few tens of nanoseconds, our measurements indicate 
that everything behaves as if this state lasted for much longer. 

To summarize our finding in few words, one might say that the data 
do not agree with the common interpretatiorj^] of the Lienard-Weichert 

2 As a reminder, the complete solution of Maxwell equations consists of linear combinations 
of retarded and advanced potentials. 
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potential for uniformly moving charges, while seem to support the idea of 
a Coulomb field carried rigidly by the electron beam. 

We would welcome any interpretation, different from the Feynman 
conjecture or the instataneous propagation, that will help understanding 
the time/space evolution of the electric field we measure. 
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